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Abstract 

Nuclear stopping in the heavy ion collisions over a beam energy range from SIS, AGS up to SPS is 
studied in the framework of the modified UrQMD transport model, in which mean field potentials of 
both formed and "pre-formed" hadrons (from string fragmentation) and medium modified nucleon- 
nucleon elastic cross sections are considered. It is found that the nuclear stopping is influenced 
by both the stiffness of the equation of state and the medium modifications of nucleon-nucleon 
cross sections at SIS energies. At the high SPS energies, the two-bump structure is shown in the 
experimental rapidity distribution of free protons, which can be understood with the consideration 
of the "pre-formed" hadron potentials. 
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I. INTRODUCTION 



Since 1980's the heavy ion collisions (HICs) in terrestrial laboratories have been becom- 
ing an important way to investigate properties of hot and dense nuclear matter []]-8]. In 

particular, the study of transport phenomena in nuclear reactions is of major importance 

n 

in the understanding of many fundamental properties [9]. And, more interest was focused 
on extracting the equation of state (EoS) of nuclear matter from the comparison of micro- 
scopic transport models with experimental measurements. Recently, the effect of medium 
modifications on two-body collisions is received more and more attention. 



In one of the attempts to obtain information about the EoS from heavy ion data [10| , it 
is made clear that progress on this topic requires improved understanding of the momentum 
dependence of mean fields generated in HICs as well as an extensive modification according 
to experimental information on the degree of stopping achieved An optimal condition 
for nuclear matter compressed to form a dense medium is that the two colliding heavy 
ions are fully stopped by each other during the process of interaction, before the system 



starts to expand 12| . Information on the stopping can be obtained by studying the rapidity 
distributions of fragments or free nucleons in the transverse and longitudinal directions. In 
the ratio of the widths of the transverse to the longitudinal rapidity distributions was 
proposed as an indicator of the stopping degree. 

The main purpose of this work is to extract the information of nuclear stopping by the 
comparison of the rapidity distributions of protons and other stopping related observable 
from a transport-model simulation with data. Meanwhile, medium modifications on inter- 
actions of particles in the dense matter can be detected as well. This goal can be achieved 
by studying the excitation function of the stopping from Au+Au collisions at SIS energies 
and the rapidity distribution of free protons from Au+Au/Pb+Pb collision at AGS and SPS 
energies, respectively, within a transport model — The Ultra-relativistic Quantum Molecular 
Dynamics (UrQMD) model. The advantages of this method are: (1) to directly compare 
existing data in each energy region, and (2) to minimize the uncertainties coming from initial 
conditions and final freeze-outs when more models are adopted. 
II. URQMD TRANSPORT MODEL 

The UrQMD model is a microscopic many-body transport approach and can be applied 



2 



to study pp, pA and AA interactions over an energy range from SIS to RHIC. This transport 
model is based on the covariant propagation of color strings, constituent quarks and diquarks 
(as string ends) accompanied by mesonic and baryonic degree of freedom 13j. In present 
model, the subhadronic degrees of freedom enter via the introduction of a formation time 
for hadrons produced in the fragmentation of strings 

HQ, which are dominant at the 



early stage of HICs with high SPS and RHIC energies. While at SIS and AGS energies, the 
new particles are produced from the decay of resonances. During the hadronic transport, it 
is known that two ingredients should be taken into account with care if a better comparison 
with data is needed: mean-field potential and two-body scattering cross section of particles 
(e.g., Ref. (rz|). 

A. The mean-field treatments 



The UrQMD model is based on parallel principles as the quantum molecular dynamics 
(QMD) model: hadrons are represented by Gaussian wave packets in phase space and the 
phase space of hadron i is propagated according to Hamilton's equation of motion 181 ] . 
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Here, r and p are the coordinate and momentum of the hadron i, respectively. The Hamil- 
tonian H consists of the kinetic energy T and the effective interaction potential energy 
U, 

H = T + U. (2) 

In the standard UrQMD model, the potential energy U includes the two-body and three- 
body Skyrme-, Yukawa-, Coulomb- and Pauli-terms 18l. Il9l|. 



i = u { Z + u^l + u Yuk + U Con + I ■ 



sky 
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For a better description of experimental data at SIS energies, more potential terms have 
to be considered [^(J. In the modified version of UrQMD (based on the version 2.0), the 
following two terms are further added: (1) the density dependent symmetry potential term 
U Bym and (2) the momentum-dependent term U m d 21[. Both the potential terms are very 
important for the dynamics of the intermediate-energy neutron-rich HICs. In this work four 
parameter sets for EoS are used for comparison: H-EoS, S-EoS, HM-EoS and SM-EoS, which 



can be found in Ref. 



20] 



At higher beam energies (AGS and SPS energies), the Yukawa-, Pauli- and symmetry- 
potentials of baryons become negligible, while the Skyrme- and the momentum-dependent 
parts of potentials still influence the whole dynamical process of HICs 22| . At SPS energies, 
the new production mechanism of particles (string fragmentation) plays more and more 
important role, in which the formation time of hadrons from the string fragmentation is 
determined by a "yo-yo" mode 18|, [l9|. During the formation time, the "pre-formed" parti- 



cles (string fragments that will be projected onto hadron states later on) are usually treated 
to be free-streaming, while reduced cross sections are only included for leading hadrons. 



In previous calculations la, |l9|, |23j, the interaction of the newly produced "pre-formed" 



particles is not taken into account. Recently, the mean-field potentials for both formed 
and "pre-formed" particles are considered for a better understanding of HBT time-related 
puzzle 24|. Meanwhile, in Ref. 24|, the rapidity distribution of net-protons from HICs at 
the SPS energy 158A GeV is shown to have a two-bump structure with the consideration 
of the "pre-formed" hadron potentials, which explains data fairly well. In this paper, more 
analyses about free protons at all SPS energies will be shown. 

At AGS and SPS energies, the relativistic effect on the relative distance and the relative 



momentum and a covariance-related reduced factor used for the update of potentials 
are considered in calculations. 
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25] 



B. The in-medium nucleon-nucleon (NN) elastic cross sections 

Besides the updates of the mean field part mentioned above, the influence of the medium 
modification on two-nucleon cross sections at the intermediate energy region should be also 
considered. In the present work we consider medium modifications on nucleon-nucleon (NN) 
elastic cross sections in the modified UrQMD model. For the inelastic channels, we still use 
the experimental free-space cross sections. It is believed that this assumption has minor 
effect on our present study at SIS energies. At present, three forms of in- medium NN elastic 
cross sections are considered, they are (1) o , the free nucleon-nucleon elastic cross section. 



(2) ai*, which is based on the extended QHD theory and reads as [26|, [27] 

( x 1 * = F( M ,a,p)a free , (4) 
where the medium correction factor F depends on the nuclear reduced density u = Pi/po, 
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the isospin-asymmetry a = (p n — p p )/pi, and the relative momentum of two colliding nuclei. 
The pi, p n and p„ are the nuclear, neutron and proton densities, respectively More explicitly, 
the factor FisflEl 



F(u,a,p) 



p p . p p 



(5) 



where 



F u p = 1 + [§ezp(-«/0.54568) - |]/[1 + fam/po)*], Pnn <1 GeV/c; 
F a p = 1 + [7^(0.85/(1 + 3.25u))a]/[l + (pnn/Po)"], Pnn <1 GeV/c; 



(6) 



- 1 a,« 



Pnn >1 GeV/c. 



Here j»nn is the relative momentum in the NN center-of-mass system; Ty = —1, +1, and in 
the case of i = j = p, i = j = n, and % ^ j, respectively; r\ is set to —1 for a nonrelativistic 
typed splitting on the proton-proton and neutron-neutron elastic cross sections in the isospin- 
asymmetric nuclear medium. The other parameters po and k, which influence the slope of the 
momentum dependence of the reduction factor F u , are still with somewhat uncertainty j^]. 
In this work, we choose po = 0.5 GeV/c and k = 6 as an example. Employing this approach, 
it was found that the in-medium NN elastic cross sections were suppressed seriously at low 
relative momenta than at higher one depending on the medium density, which is similar to 



the Brueckner relativistic approach 28 



29[. (3) a"2*, as in Ref. 30], which reads as 



&2 



(1 - ^V free , (7) 

where £ = 0.5 for < 0.25A GeV in this work. It is easy to find that the momentum 
constraint is not considered in 02*. Further, the density dependence of is stronger than 
that of <7i*. 



For calculations at SIS, a conventional phase-space coalescence model [31] is used to 
construct clusters, in which nucleons with relative distances smaller than R and relative 
momenta smaller than P Q are considered to belong to one cluster. Fig. [TJ shows normalized 
rapidity distributions of fragments with proton number Z=l, 3, and 8 (from top to bottom 
panels) in the longitudinal (left panel) and transverse (right panel) directions for central 
Au+Au collisions at 0.15^4 GeV. Two (Rq, Pq) parameter sets, (3.5 fm, 0.2 GeV/c) and (3.0 
fm, 0.2 GeV/c), are adopted in the calculations. The results are shown with lines and the 



FOPI data 



12] are shown by scattered stars. It seems that the parameter set (3.0 fm, 0.2 



GeV/c) gives a better description of the FOPI data. Therefore, this parameter set is used 



Au+Au, E, =0.1 5A GeV, b=0-2 fm 




FIG. 1: Normalized rapidity distributions of fragments with proton number Z=l, 3, and 8 (from 
top to bottom panels) in the longitudinal (left panel) and transverse (right panel) directions for 
central Au+Au collisions at 0.15^4 GeV. Two (Rq, Pq) parameter sets, (3.5 fm, 0.2 GeV/c) and 
(3.0 fm, 0.2 GeV/c), are adopted in the coalescence- model calculations, which are shown with lines. 
FOPI data [3] are shown by scattered stars. 

in the following calculations at SIS energies in this work. While at AGS and SPS energies, 
the coalescence model is not used as usual (partly because of the rich production of new 
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baryons) so that all nucleons at freeze-out are taken to be free. 

III. NUCLEAR STOPPING AND THE RAPIDITY DISTRIBUTIONS 



A. vartl at SIS energies 

As a measure of the nuclear stopping degree {fl, the FOPI Collaboration introduced 
a new observable vartl which was defined by the ratio of the variances of the transverse 
to the longitudinal rapidity distributions of fragments. For central Au+Au collisions, it is 
found that the rapidity distributions in the x and y directions are nearly the same, thus the 
transverse rapidity distributions are plotted approximately with the rapidity distributions 
in the x direction. Numerically, the vartl is defined as 

vartl = T dN / d y x /T dN / dyz , (8) 

where T dN / dyx (T dN / dyz ) is the width of the rapidity distribution of fragments in the x (z) 
direction and reads as 



F d N/dy x>z = y(yl, z ), (9) 

/ 2 \ (yx,zNy x z ) , , 

(y x ,z) = 77 • (io) 

Here N Vx z and N a \\ are yields of fragments in each y x (or y z ) rapidity bin and in the whole 
rapidity region, respectively. It is easy to understand that vartl < 1 stands for an incom- 
plete stopping or nuclear transparency, and vartl > 1 for a strong transverse expansion or 
collectivity. Obviously, vartl = 1 when a full stopping occurs. 

The excitation function of vartl for central Au+Au collisions is shown in Fig. [2] within 



the beam energy region 0.09v4 — 1.5A GeV . The FOPI data rjj are shown by stars while 
the UrQMD calculations are shown by lines with symbols. The vartl value is calculated for 
fragments with the proton number Z < 10. In the calculations, results with the cascade 
mode and with various EoS are shown. The free NN cross sections are adopted in the 
calculations. It is seen that the vartl value of the cascade mode is always less than 1 
and decreases monotonously with the increase of beam energy, which implies less and less 
stopping strength in the system. At ~ 0.3 A — 1A GeV calculated values of vartl are 
smaller than data while it is larger than data at lower beam energies. When the mean field 
is considered, the potentials reinforce the bound of nucleons and a stronger collectivity is 
shown in the transverse direction. Among the calculations with EoS, softer EoS gives a 
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+IG. 2: Excitation function of vartl for central Au+Au collisions at SIS energies. The FOPI data 
111 ] are shown by stars while the UrQMD calculations with various EoS are shown by lines with 
symbols. 

smaller vartl value while the momentum dependent term in the potential plays a negligible 
role. We also find that only a soft EoS can not describe the excitation function of the FOPI 
data without considering medium modifications of two-body collisions. Next, based on the 
result with the SM-EoS, we will give a further investigation of the effect of the medium 
modifications of NN elastic cross sections on the vartl. 

Fig.|3]iilustrates the calculated excitation function of vartl with the medium modified NN 
elastic cross section a\ as well as the free one <r free . It is seen clearly that a large reduction 
of cross sections at lower beam energies leads to obvious transparency so that the calculated 
vartl with a\ are largely decreased at low SIS energies. While at high SIS energies the vartl 
value is much less affected and slightly higher than data. As mentioned in Eq. [6j this might 
be due to the fixed pnn cut adopted. We would not modify this just for fitting data since 
the medium modifications on inelastic channels are still an open question. We just wish to 
stress the importance of medium modifications of cross sections on the nuclear stopping at 
moderate SIS energies. 

For i?i a b < 0.25A GeV, the results with <j\ are still higher than data which implies that 
a stronger reduction factor on the elastic cross sections is required. Fig. H] further shows the 
calculation with 02* (with a stronger reduction factor on the NN elastic cross section, as 
seen in Eq. [7j) for < 0.25A GeV. The comparison with data is fairly well and same as 
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FIG. 3: Excitation function of fart/ with the medium modified NN elastic cross section a\* as 
well as the free one <r free . The SM-EoS is adopted in calculations. FOPI data 
comparison. 



Ill ] are shown for 
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done in Ref. [17]. 

B. Rapidity distribution at AGS and SPS energies 



At AGS and SPS energies, as the rapidity distribution of fragments in the transverse 
direction has not been provided by experiments, we study the nuclear stopping with the 
longitudinal rapidity distribution. Figs. [5] and [6] depict the rapidity distributions of protons 
for central Au+Au collisions at AGS and for central Pb+Pb collisions at SPS (< 5% of 
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total cross section <tt), respectively. The (preliminary) experimental data of free protons 



are taken from 



321-1351] . In the calculations, besides a cascade mode shown in the left panel, 



we also show the results with potentials of both formed and "pre-formed" hadrons ( "pf-part 
& f-B SM-EoS" ) in the right panel. Cross sections used in the model are not modified by the 
nuclear medium in this energy region. Since protons belonging to fragments are included 
in calculations of the rapidity distribution, the calculation results of the proton number 



are somewhat larger than data, especia. 
[6] as well as in previous calculations 13 



ly at low beam energies as shown in Figs. [5] and 



36|. In Fig. [5j it is found that the shape of the 



rapidity distributions of measured protons changes from one peak at mid-rapidity with no 
shoulder to two shoulders when increasing beam energy from 2A GeV to 11A GeV. The 
cascade calculations always give a Gaussian-like distribution at y < 1.0, while calculations 
with potentials are much closer to data. With the increase of beam energy from AGS to 
SPS, the experimental rapidity distribution changes further to a plateau and finally to a 
two-bump structure. Again, the calculations with cascade mode cannot describe the shape 
of the rapidity distribution of protons completely. The stronger repulsion at early stage 
introduced by potentials makes a wider rapidity distribution of protons in the longitudinal 



direction 



241 ] . The gap of two peaks becomes wider with the increase of beam energy 



Especially, at 160A GeV the rapidity distribution of protons shows clearly two peaks at 
y ~ 1.5. These features can be reasonably reproduced by the calculations with both the 
formed and "pre-formed" hadron potentials shown in the right panel of Fig. |6j 

We also calculate the rapidity distribution of emitted As for central Pb+Pb collisions at 
AOA GeV and 160A GeV with and without formed and "pre-formed" hadron potentials as 
shown in Fig. [7J Calculations with and without potentials (lines) are compared to the NA49 
data |37| (stars). Same as data, the yields represent the sum A + S°. It is seen clearly that 
calculations with potentials are in good agreement with data at both beam energies, which 
is due to a larger transparency introduced by the strongly repulsive mean field at the early 
stage. As is known that at the AGS and SPS energies the yields of hyperons are somewhat 



overestimated in the UrQMD cascade calculations with version less than 2.1 jlSL Il9l |23| . 
which is also shown in Fig. [7J In order to solve this problem, alternatively, starting from 
the version 2.1 (and the recently published v2.3), the UrQMD group considers additional 
high mass resonances that are explicitly produced and propagated in s-channel processes 



with invariant masses up to yfs < 3 GeV[23j, [38(. This treatment leads to lower yield of the 
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FIG. 5: Rapidity distributions of protons at AGS energies 2A, AA, 6A, 8A, and 11A GeV for 
central Au+Au collisions. Calculations with cascade (left panel) and with potentials "pf-part & 
f-B SM-EoS" (right panel) are shown with lines. Experimental data of free protons taken from 



E895 



32fl and E802 



331 ] Collaborations are shown with scattered symbols. 



strange particles so that a nice agreement with A data from central Pb+Pb collisions at SPS 



energies was also shown in previous calculations [131. 138l|. Therefore, it deserves much more 
investigations to deeply understand the effects of mean field potentials and the decay of high 
mass resonances on, e.g., particle production and collective flows, which are in progress. 
IV. SUMMARY AND OUTLOOK 



In summary, we have presented the excitation function of the nuclear stopping described 
by vartl of light fragments for central Au+Au reactions with beam energies from 0.09/1 GeV 
to 1.5A GeV and the rapidity distribution of protons and As for central Au+Au/Pb+Pb 
reactions in the energy region 2 A — 160 A GeV. The modified UrQMD transport model (based 
on the version 2.0) has been used in all calculations. Based on the model we investigate the 
effects of both the mean-field potentials and medium modifications of nucleon-nucleon elastic 
cross sections on the nuclear stopping under the same initial and final freeze-out conditions. 
It is found that the nuclear stopping is influenced by both the stiffness of the equation of 
state and the medium modifications of nucleon-nucleon elastic cross sections for reactions 
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FIG. 6: Rapidity distributions of protons at SPS energies 20A, 30A, 40A, 80A, and 160 A GeV for 
central Pb+Pb collisions. Calculations with cascade (left panel) and with potentials "pf-part & 
-B SM-EoS" (right panel) are shown by lines. Preliminary data of free protons taken from NA49 
34 . 3^ ] Collaboration are shown by scattered symbols. 




FIG. 7: Rapidity distributions of As from central Pb+Pb collisions at 40A GeV (left panel) and 
160^4 GeV (right panel). Calculations with and without potentials (lines) are compared to the 
NA49 data [37] (stars). The open stars are data points reflected around mid-rapidity. 
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at SIS energies. And it reaches a well defined plateau of maximal stopping centered around 
(0.5 ± 0.3) A GeV with a fast drop on both sides. At AGS and SPS energies, the degree of 
nuclear stopping decreases continuously. In the high SPS energy region, as the transparency 
of matter is high, the two-bump structure is shown in the experimental rapidity distribution 
of free protons in the longitudinal direction. Our calculations show that considering the 
potentials of both formed and "pre-formed" hadrons can improve the agreement between 
calculation results and data. But the form of the potentials is still simple and rough and 
further improvement is needed. The work on this aspect is underway. 
Acknowledgements 

We acknowledge support by the Frankfurt Center for Scientific Computing (CSC). This 
work is supported by: the Key Project of the Ministry of Education of China under Grant 
No. 209053, the National Natural Science Foundation of China under Grant Nos. 10675077, 
10975095, 10675172, 10875031, 10905021, 10979023, the National Basic Research Program 
of China under Grant No. 2007CB209900, and the Natural Science Foundation of Zhejiang 
Province under grant No. Y6090210. 



13 



[1] J. Randrup, Nucl. Phys. A 314, 429 (1979). 
[2] Z. Li, et al, Nucl. Phys. A 559, 603 (1993). 
[3] Z. Li, et al, J. Phys. G 20, 1829 (1994). 

[4] T. K. Choi, M. Maruyama and F. Takagi, Phys. Rev. C 55, 848 (1997). 

[5] J. D. Bjorken, Phys. Rev. D 27, 140 (1983). 

[6] Q. Li and Z. Li, Chin. Phys. Lett. 19, 321 (2002). 

[7] H. Petersen, Q. Li, X. Zhu and M. Bleicher, Phys. Rev. C 74, 064908 (2006). 

[8] P. Danielewicz, B. Barker and L. Shi, AIP Conf. Proc. 1128, 104 (2009). 

[9] C. Escano- Rodriguez et al. [INDRA Collaboration and ALADIN Collaboration], 
arXiv:nucl-ex/0503007 
[10] P. Danielewicz, R. Lacey and W. G. Lynch, Science 298, 1592 (2002). 
[11] W. Reisdorf et al. [FOPI Collaboration], Phys. Rev. Lett. 92, 232301 (2004). 
[12] A. Andronic, J. Lukasik, W. Reisdorf and W. Trautmann, Eur. Phys. J. A 30, 31 (2006). 
[13] H. Petersen, M. Bleicher, S. A. Bass and H. Stocker, arXiv:0805.0567 [hep-ph]. 
[14] B. Andersson, G. Gustafson and B. Nilsson-Almqvist, Nucl. Phys. B 281, 289 (1987). 
[15] B. Nilsson-Almqvist and E. Stenlund, Comput. Phys. Commun. 43, 387 (1987). 
[16] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 
[17] Y. Zhang, Z. Li and P. Danielewicz, Phys. Rev. C 75, 034615 (2007). 
[18] S. A. Bass et al, Prog. Part. Nucl. Phys. 41, 255 (1998). 
[19] M. Bleicher et al, J. Phys. G 25, 1859 (1999). 

[20] Q. Li, Z. Li, S. Soff, M. Bleicher and H. Stoecker, J. Phys. G 32, 151 (2006). 

[21] S. A. Bass, C. Hartnack, H. Stoecker and W. Greiner, Phys. Rev. C 51, 3343 (1995). 

[22] Q. Li and M. Bleicher, J. Phys. G 36, 015111 (2009). 

[23] E. L. Bratkovskaya et al, Phys. Rev. C 69, 054907 (2004). 

[24] Q. Li, M. Bleicher and H. Stocker, Phys. Lett. B 659, 525 (2008). 

[25] M. Isse, A. Ohnishi, N. Otuka, P. K. Sahu and Y. Nara, Phys. Rev. C 72, 064908 (2005). 

[26] Q. Li, Z. Li, S. Soff, M. Bleicher and H. Stoecker, J. Phys. G 32, 407 (2006). 

[27] Q. Li, C. Shen and M. Di Toro. larXiv:0908.2825 [nucl-th]. 

[28] C. Fuchs, A. Faessler and M. El-Shabshiry, Phys. Rev. C 64, 024003 (2001). 



14 



[29] T. Gaitanos, C. Puchs and H. H. Wolter, Phys. Lett. B 609, 241 (2005). 
[30] D. Klakow, G. Welke and W. Bauer, Phys. Rev. C 48, 1982 (1993). 

[31] H. Kruse, B. V. Jacak, J. J. Molitoris, G. D. Westfall and H. Stoecker, Phys. Rev. C 31, 1770 
(1985). 

[32] J. L. Klay et al. [E895 Collaboration], Phys. Rev. Lett. 88, 102301 (2002). 
[33] Y. Akiba et al. [E802 Collaboration], Nucl. Phys. A 610 (1996) 139C. 

[34] C. Blume [NA49 Collaboration], J. Phys. G 34, S951 (2007); private communication of pre- 
liminary data. 

[35] H. Strobele [NA49 Collaboration], |ar2avi0908 i 2777J [nucl-ex]. 

[36] W. Reisdorf, Prog. Theor. Phys. Suppl. 140, 111 (2000). 

[37] T. Anticic et al. [NA49 Collaboration], Phys. Rev. C 80, 034906 (2009). 

[38] H. Petersen, M. Mitrovski, T. Schuster and M. Bleicher, Phys. Rev. C 80, 054910 (2009). 



15 



Transport model study of nuclear stopping in heavy ion collisions 
over an energy range from 0.09^4 GeV to 160A GeV 

Ying Yuan, 1,3 * Qingfeng Li, 2 ^ Zhuxia Li, 3 * an d Fu-Hu Liu 1§ 

1 ) Institute of Theoretical Physics, Shanxi University, Taiyuan, Shanxi 030006, China 
2) School of Science, Huzhou Teachers College, Huzhou, Zhejiang 313000, China 

3) China Institute of Atomic Energy, 
P. O. Box 275 (18), Beijing 102413, China 

Abstract 

Nuclear stopping in the heavy ion collisions over a beam energy range from SIS, AGS up to SPS is 
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by both the stiffness of the equation of state and the medium modifications of nucleon-nucleon 
cross sections at SIS energies. At the high SPS energies, the two-bump structure is shown in the 
experimental rapidity distribution of free protons, which can be understood with the consideration 
of the "pre-formed" hadron potentials. 
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I. INTRODUCTION 



Since 1980's the heavy ion collisions (HICs) in terrestrial laboratories have been becom- 
ing an important way to investigate properties of hot and dense nuclear matter []]-8]. In 

particular, the study of transport phenomena in nuclear reactions is of major importance 

n 

in the understanding of many fundamental properties [9]. And, more interest was focused 
on extracting the equation of state (EoS) of nuclear matter from the comparison of micro- 
scopic transport models with experimental measurements. Recently, the effect of medium 
modifications on two-body collisions is received more and more attention. 



In one of the attempts to obtain information about the EoS from heavy ion data [10| , it 
is made clear that progress on this topic requires improved understanding of the momentum 
dependence of mean fields generated in HICs as well as an extensive modification according 
to experimental information on the degree of stopping achieved An optimal condition 
for nuclear matter compressed to form a dense medium is that the two colliding heavy 
ions are fully stopped by each other during the process of interaction, before the system 



starts to expand 12| . Information on the stopping can be obtained by studying the rapidity 
distributions of fragments or free nucleons in the transverse and longitudinal directions. In 
the ratio of the widths of the transverse to the longitudinal rapidity distributions was 
proposed as an indicator of the stopping degree. 

The main purpose of this work is to extract the information of nuclear stopping by the 
comparison of the rapidity distributions of protons and other stopping related observable 
from a transport-model simulation with data. Meanwhile, medium modifications on inter- 
actions of particles in the dense matter can be detected as well. This goal can be achieved 
by studying the excitation function of the stopping from Au+Au collisions at SIS energies 
and the rapidity distribution of free protons from Au+Au/Pb+Pb collision at AGS and SPS 
energies, respectively, within a transport model — The Ultra-relativistic Quantum Molecular 
Dynamics (UrQMD) model. The advantages of this method are: (1) to directly compare 
existing data in each energy region, and (2) to minimize the uncertainties coming from initial 
conditions and final freeze-outs when more models are adopted. 
II. URQMD TRANSPORT MODEL 

The UrQMD model is a microscopic many-body transport approach and can be applied 



2 



to study pp, pA and AA interactions over an energy range from SIS to RHIC. This transport 
model is based on the covariant propagation of color strings, constituent quarks and diquarks 
(as string ends) accompanied by mesonic and baryonic degree of freedom 13j. In present 
model, the subhadronic degrees of freedom enter via the introduction of a formation time 
for hadrons produced in the fragmentation of strings 

HQ, which are dominant at the 



early stage of HICs with high SPS and RHIC energies. While at SIS and AGS energies, the 
new particles are produced from the decay of resonances. During the hadronic transport, it 
is known that two ingredients should be taken into account with care if a better comparison 
with data is needed: mean-field potential and two-body scattering cross section of particles 
(e.g., Ref. (rz|). 

A. The mean-field treatments 



The UrQMD model is based on parallel principles as the quantum molecular dynamics 
(QMD) model: hadrons are represented by Gaussian wave packets in phase space and the 
phase space of hadron i is propagated according to Hamilton's equation of motion 181 ] . 

■ dH - OH 



Pi 



(1) 



Here, r and p are the coordinate and momentum of the hadron i, respectively. The Hamil- 
tonian H consists of the kinetic energy T and the effective interaction potential energy 
U, 

H = T + U. (2) 

In the standard UrQMD model, the potential energy U includes the two-body and three- 
body Skyrme-, Yukawa-, Coulomb- and Pauli-terms 18l. Il9l|. 



i = u { Z + u^l + u Yuk + U Con + I ■ 



sky 



pau' 



(3) 



For a better description of experimental data at SIS energies, more potential terms have 
to be considered [^(J. In the modified version of UrQMD (based on the version 2.0), the 
following two terms are further added: (1) the density dependent symmetry potential term 
U Bym and (2) the momentum-dependent term U m d 21[. Both the potential terms are very 
important for the dynamics of the intermediate-energy neutron-rich HICs. In this work four 
parameter sets for EoS are used for comparison: H-EoS, S-EoS, HM-EoS and SM-EoS, which 



can be found in Ref. 



20] 



At higher beam energies (AGS and SPS energies), the Yukawa-, Pauli- and symmetry- 
potentials of baryons become negligible, while the Skyrme- and the momentum-dependent 
parts of potentials still influence the whole dynamical process of HICs 22| . At SPS energies, 
the new production mechanism of particles (string fragmentation) plays more and more 
important role, in which the formation time of hadrons from the string fragmentation is 
determined by a "yo-yo" mode 18|, [l9|. During the formation time, the "pre-formed" parti- 



cles (string fragments that will be projected onto hadron states later on) are usually treated 
to be free-streaming, while reduced cross sections are only included for leading hadrons. 



In previous calculations la, |l9|, |23j, the interaction of the newly produced "pre-formed" 



particles is not taken into account. Recently, the mean-field potentials for both formed 
and "pre-formed" particles are considered for a better understanding of HBT time-related 
puzzle 24|. Meanwhile, in Ref. 24|, the rapidity distribution of net-protons from HICs at 
the SPS energy 158A GeV is shown to have a two-bump structure with the consideration 
of the "pre-formed" hadron potentials, which explains data fairly well. In this paper, more 
analyses about free protons at all SPS energies will be shown. 

At AGS and SPS energies, the relativistic effect on the relative distance and the relative 



momentum and a covariance-related reduced factor used for the update of potentials 
are considered in calculations. 



22 



25] 



B. The in-medium nucleon-nucleon (NN) elastic cross sections 

Besides the updates of the mean field part mentioned above, the influence of the medium 
modification on two-nucleon cross sections at the intermediate energy region should be also 
considered. In the present work we consider medium modifications on nucleon-nucleon (NN) 
elastic cross sections in the modified UrQMD model. For the inelastic channels, we still use 
the experimental free-space cross sections. It is believed that this assumption has minor 
effect on our present study at SIS energies. At present, three forms of in- medium NN elastic 
cross sections are considered, they are (1) o , the free nucleon-nucleon elastic cross section. 



(2) ai*, which is based on the extended QHD theory and reads as [26|, [27] 

( x 1 * = F( M ,a,p)a free , (4) 
where the medium correction factor F depends on the nuclear reduced density u = Pi/po, 
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the isospin-asymmetry a = (p n — p p )/pi, and the relative momentum of two colliding nuclei. 
The pi, p n and p„ are the nuclear, neutron and proton densities, respectively More explicitly, 
the factor FisflEl 



F(u,a,p) 



p p . p p 



(5) 



where 



F u p = 1 + [§ezp(-«/0.54568) - |]/[1 + fam/po)*], Pnn <1 GeV/c; 
F a p = 1 + [7^(0.85/(1 + 3.25u))a]/[l + (pnn/Po)"], Pnn <1 GeV/c; 



(6) 



- 1 a,« 



Pnn >1 GeV/c. 



Here j»nn is the relative momentum in the NN center-of-mass system; Ty = —1, +1, and in 
the case of i = j = p, i = j = n, and % ^ j, respectively; r\ is set to —1 for a nonrelativistic 
typed splitting on the proton-proton and neutron-neutron elastic cross sections in the isospin- 
asymmetric nuclear medium. The other parameters po and k, which influence the slope of the 
momentum dependence of the reduction factor F u , are still with somewhat uncertainty j^]. 
In this work, we choose po = 0.5 GeV/c and k = 6 as an example. Employing this approach, 
it was found that the in-medium NN elastic cross sections were suppressed seriously at low 
relative momenta than at higher one depending on the medium density, which is similar to 



the Brueckner relativistic approach 28 



29[. (3) a"2*, as in Ref. 30], which reads as 



&2 



(1 - ^V free , (7) 

where £ = 0.5 for < 0.25A GeV in this work. It is easy to find that the momentum 
constraint is not considered in 02*. Further, the density dependence of is stronger than 
that of <7i*. 



For calculations at SIS, a conventional phase-space coalescence model [31] is used to 
construct clusters, in which nucleons with relative distances smaller than R and relative 
momenta smaller than P Q are considered to belong to one cluster. Fig. [TJ shows normalized 
rapidity distributions of fragments with proton number Z=l, 3, and 8 (from top to bottom 
panels) in the longitudinal (left panel) and transverse (right panel) directions for central 
Au+Au collisions at 0.15^4 GeV. Two (Rq, Pq) parameter sets, (3.5 fm, 0.2 GeV/c) and (3.0 
fm, 0.2 GeV/c), are adopted in the calculations. The results are shown with lines and the 



FOPI data 



12] are shown by scattered stars. It seems that the parameter set (3.0 fm, 0.2 



GeV/c) gives a better description of the FOPI data. Therefore, this parameter set is used 



Au+Au, E, =0.15/1 GeV, b=0-2 fm 




FIG. 1: Normalized rapidity distributions of fragments with proton number Z=l, 3, and 8 (from 
top to bottom panels) in the longitudinal (left panel) and transverse (right panel) directions for 
central Au+Au collisions at 0.15^4 GeV. Two (Rq, Pq) parameter sets, (3.5 fm, 0.2 GeV/c) and 
(3.0 fm, 0.2 GeV/c), are adopted in the coalescence- model calculations, which are shown with lines. 
FOPI data [3] are shown by scattered stars. 

in the following calculations at SIS energies in this work. While at AGS and SPS energies, 
the coalescence model is not used as usual (partly because of the rich production of new 
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baryons) so that all nucleons at freeze-out are taken to be free. 

III. NUCLEAR STOPPING AND THE RAPIDITY DISTRIBUTIONS 



A. vartl at SIS energies 

As a measure of the nuclear stopping degree {fl, the FOPI Collaboration introduced 
a new observable vartl which was defined by the ratio of the variances of the transverse 
to the longitudinal rapidity distributions of fragments. For central Au+Au collisions, it is 
found that the rapidity distributions in the x and y directions are nearly the same, thus the 
transverse rapidity distributions are plotted approximately with the rapidity distributions 
in the x direction. Numerically, the vartl is defined as 

vartl = T dN / d y x /T dN / dyz , (8) 

where T dN / dyx (T dN / dyz ) is the width of the rapidity distribution of fragments in the x (z) 
direction and reads as 



F d N/dy x>z = y(yl, z ), (9) 

/ 2 \ (yx,zNy x z ) , , 

(y x ,z) = 77 • (io) 

Here N Vx z and N a \\ are yields of fragments in each y x (or y z ) rapidity bin and in the whole 
rapidity region, respectively. It is easy to understand that vartl < 1 stands for an incom- 
plete stopping or nuclear transparency, and vartl > 1 for a strong transverse expansion or 
collectivity. Obviously, vartl = 1 when a full stopping occurs. 

The excitation function of vartl for central Au+Au collisions is shown in Fig. [2] within 



the beam energy region 0.09v4 — 1.5A GeV . The FOPI data rjj are shown by stars while 
the UrQMD calculations are shown by lines with symbols. The vartl value is calculated for 
fragments with the proton number Z < 10. In the calculations, results with the cascade 
mode and with various EoS are shown. The free NN cross sections are adopted in the 
calculations. It is seen that the vartl value of the cascade mode is always less than 1 
and decreases monotonously with the increase of beam energy, which implies less and less 
stopping strength in the system. At ~ 0.3 A — 1A GeV calculated values of vartl are 
smaller than data while it is larger than data at lower beam energies. When the mean field 
is considered, the potentials reinforce the bound of nucleons and a stronger collectivity is 
shown in the transverse direction. Among the calculations with EoS, softer EoS gives a 
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+IG. 2: Excitation function of vartl for central Au+Au collisions at SIS energies. The FOPI data 
111 ] are shown by stars while the UrQMD calculations with various EoS are shown by lines with 
symbols. 

smaller vartl value while the momentum dependent term in the potential plays a negligible 
role. We also find that only a soft EoS can not describe the excitation function of the FOPI 
data without considering medium modifications of two-body collisions. Next, based on the 
result with the SM-EoS, we will give a further investigation of the effect of the medium 
modifications of NN elastic cross sections on the vartl. 

Fig.|3]iilustrates the calculated excitation function of vartl with the medium modified NN 
elastic cross section a\ as well as the free one <r free . It is seen clearly that a large reduction 
of cross sections at lower beam energies leads to obvious transparency so that the calculated 
vartl with a\ are largely decreased at low SIS energies. While at high SIS energies the vartl 
value is much less affected and slightly higher than data. As mentioned in Eq. [6j this might 
be due to the fixed pnn cut adopted. We would not modify this just for fitting data since 
the medium modifications on inelastic channels are still an open question. We just wish to 
stress the importance of medium modifications of cross sections on the nuclear stopping at 
moderate SIS energies. 

For i?i a b < 0.25A GeV, the results with <j\ are still higher than data which implies that 
a stronger reduction factor on the elastic cross sections is required. Fig. H] further shows the 
calculation with 02* (with a stronger reduction factor on the NN elastic cross section, as 
seen in Eq. [7j) for < 0.25A GeV. The comparison with data is fairly well and same as 
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FIG. 3: Excitation function of fart/ with the medium modified NN elastic cross section a\* as 
well as the free one <r free . The SM-EoS is adopted in calculations. FOPI data 
comparison. 



Ill ] are shown for 
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SM-EoS is adopted in calculations. 



done in Ref. [17]. 

B. Rapidity distribution at AGS and SPS energies 



At AGS and SPS energies, as the rapidity distribution of fragments in the transverse 
direction has not been provided by experiments, we study the nuclear stopping with the 
longitudinal rapidity distribution. Figs. [5] and [6] depict the rapidity distributions of protons 
for central Au+Au collisions at AGS and for central Pb+Pb collisions at SPS (< 5% of 
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total cross section <tt), respectively. The (preliminary) experimental data of free protons 



are taken from 



321-1351] . In the calculations, besides a cascade mode shown in the left panel, 



we also show the results with potentials of both formed and "pre-formed" hadrons ( "pf-part 
& f-B SM-EoS" ) in the right panel. Cross sections used in the model are not modified by the 
nuclear medium in this energy region. Since protons belonging to fragments are included 
in calculations of the rapidity distribution, the calculation results of the proton number 



are somewhat larger than data, especia. 
[6] as well as in previous calculations 13 



ly at low beam energies as shown in Figs. [5] and 



36|. In Fig. [5j it is found that the shape of the 



rapidity distributions of measured protons changes from one peak at mid-rapidity with no 
shoulder to two shoulders when increasing beam energy from 2A GeV to 11A GeV. The 
cascade calculations always give a Gaussian-like distribution at y < 1.0, while calculations 
with potentials are much closer to data. With the increase of beam energy from AGS to 
SPS, the experimental rapidity distribution changes further to a plateau and finally to a 
two-bump structure. Again, the calculations with cascade mode cannot describe the shape 
of the rapidity distribution of protons completely. The stronger repulsion at early stage 
introduced by potentials makes a wider rapidity distribution of protons in the longitudinal 



direction 



241 ] . The gap of two peaks becomes wider with the increase of beam energy 



Especially, at 160A GeV the rapidity distribution of protons shows clearly two peaks at 
y ~ 1.5. These features can be reasonably reproduced by the calculations with both the 
formed and "pre-formed" hadron potentials shown in the right panel of Fig. |6j 

We also calculate the rapidity distribution of emitted As for central Pb+Pb collisions at 
AOA GeV and 160A GeV with and without formed and "pre-formed" hadron potentials as 
shown in Fig. [7J Calculations with and without potentials (lines) are compared to the NA49 
data |37| (stars). Same as data, the yields represent the sum A + S°. It is seen clearly that 
calculations with potentials are in good agreement with data at both beam energies, which 
is due to a larger transparency introduced by the strongly repulsive mean field at the early 
stage. As is known that at the AGS and SPS energies the yields of hyperons are somewhat 



overestimated in the UrQMD cascade calculations with version less than 2.1 jlSL Il9l |23| . 
which is also shown in Fig. [7J In order to solve this problem, alternatively, starting from 
the version 2.1 (and the recently published v2.3), the UrQMD group considers additional 
high mass resonances that are explicitly produced and propagated in s-channel processes 



with invariant masses up to yfs < 3 GeV[23j, [38(. This treatment leads to lower yield of the 



10 




FIG. 5: Rapidity distributions of protons at AGS energies 2A, AA, 6A, 8A, and 11A GeV for 
central Au+Au collisions. Calculations with cascade (left panel) and with potentials "pf-part & 
f-B SM-EoS" (right panel) are shown with lines. Experimental data of free protons taken from 



E895 32 and E802 



331 ] Collaborations are shown with scattered symbols. 



strange particles so that a nice agreement with A data from central Pb+Pb collisions at SPS 
energies was also shown in previous calculations [3, 



. Therefore, it deserves much more 



investigations to deeply understand the effects of mean field potentials and the decay of high 
mass resonances on, e.g., particle production and collective flows, which are in progress. 
IV. SUMMARY AND OUTLOOK 

In summary, we have presented the excitation function of the nuclear stopping described 
by vartl of light fragments for central Au+Au reactions with beam energies from 0.09/1 GeV 
to 1.5A GeV and the rapidity distribution of protons and As for central Au+Au/Pb+Pb 
reactions in the energy region 2 A — 160 A GeV. The modified UrQMD transport model (based 
on the version 2.0) has been used in all calculations. Based on the model we investigate the 
effects of both the mean-field potentials and medium modifications of nucleon-nucleon elastic 
cross sections on the nuclear stopping under the same initial and final freeze-out conditions. 
It is found that the nuclear stopping is influenced by both the stiffness of the equation of 
state and the medium modifications of nucleon-nucleon elastic cross sections for reactions 
at SIS energies. And it reaches a well defined plateau of maximal stopping centered around 
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FIG. 6: Rapidity distributions of protons at SPS energies 20A, 30A, 40A, 80A, and 160A GeV for 
central Pb+Pb collisions. Calculations with cascade (left panel) and with potentials "pf-part & 
:'-B SM-EoS" (right panel) are shown by lines. Preliminary data of free protons taken from NA49 
34 . Collaboration are shown by scattered symbols. 



(0.5 ± 0.3)A GeV with a fast drop on both sides. At AGS and SPS energies, the degree of 
nuclear stopping decreases continuously. In the high SPS energy region, as the transparency 
of matter is high, the two-bump structure is shown in the experimental rapidity distribution 
of free protons in the longitudinal direction. Our calculations show that considering the 
potentials of both formed and "pre-formed" hadrons can improve the agreement between 
calculation results and data. But the form of the potentials is still simple and rough and 
further improvement is needed. The work on this aspect is underway. 
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